Non-adiabatic Arbitrary Geometric Phase Gate in 2-qubit NMR Model 
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We study a 2-qubit nuclear spin system for realizing an arbitrary geometric quantum phase gate 
by means of non-adiabatic operation. A single magnetic pulse with multi harmonic frequencies is 
applied to manipulate the quantum states of 2-qubit instantly. Using resonant transition approxi- 
mation, the time dependent Hamiltonian of two nuclear spins can be solved analytically. The time 
evolution of the wave function is obtained without adiabatic approximation. The parameters of 
magnetic pulse, such as the frequency, amplitude, phase of each harmonic part as well as the time 
duration of the pulse, are determined for achieving an arbitrary non-adiabatic geometric phase gate. 
The derivation of non-adiabatic geometric controlled phase gates and A-A phase are also addressed. 



Taking the advantage of quantum superposition and 
entanglement, a powerful new computational algorithms, 
such as factoring large number— and searching an un- 
sorted database 2 -, has been created. Further study of 
quantum information processing and realization of quan- 
tum computer has attracted large numbers of theoret- 
ical and experimental scientists. For implementation 
of quantum computing, the ability to perform any sin- 
gle quantum bit (qubit) rotations (SU(2) transforma- 
tion) and 2-qubit controlled operations, i.e. controlled 
not (CNOT) or controlled phase gate (CPG), is the key 
requiremenl£ii£i£i£i£i2. Moreover, each individual qubit 
should have a long decoherence time to prevent the de- 
formation of quantum information in qubits for enough 
long time. In order to implement a controlled quantum 
gate (CNOT or CPG), the "interaction" between qubits 
is necessary. In this letter, we are interested in the phase 
gate of two nuclear spin system. In general, the phase 
shift of the quantum state stored in qubits after a cyclic 
time evolution includes both the dynamical and geomet- 
ric part. In some condition, the phase can be dependent 
only on the path of the state evolution, not the dynam- 
ics at all. In the case, the relevant transformation of the 
quantum state is called as a geometric phase gat o 10 ' n i 12 
which has been discussed and implemented in many 
schemes including trapped ion o 13 ' 14 , cavity QE D 15 ' 16 i 17 
and nuclear magnetic resonanc e) 18 : 19 ' 20 . However, the 
adiabatic implementation of geometric phase gate may 
create some considerable errors^ 1 -. Non-adiabatic imple- 
mentation of quantum information procession is required 
in the practice since any quantum computing is requested 
as fast as possible. Hence, the theoretical study of non- 
adiabatic process in quantum information procession is 



important and has attracted many scientists. In this pa- 
per, we focus on the study of the implementation of non- 
adiabatic geometric phase gate 2 ^. 

Since the nuclear spins are well isolated from the en- 
vironment, their decoherence time are far longer. Mean- 
while, nuclei with spin 1/2 are the natural as qubits in 
quantum information processing. In terms of nuclear 
magnetic resonance (NMR) , the system of nuclear spins 
has been a good candidate to demonstrate the quantum 
computational algorithm a 23 ' 24 ! 25 . Over the recent several 
years, many kinds of complex quantum information pro- 
cessing have been realized by using NMR, ranging from 
two to seven qubits in size, in liquid sample o 26 ' 27 : 28 ' 29 : 30 
and in solid state sample o 31 i 32 ' 33 . The physical manipu- 
lation of nuclear spin state can be realized by sequences 
of magnetic pulses with some active resonant frequencies. 

In this paper, we study the 2-qubit system that is com- 
posed by two weak coupled nuclear spins in an external 
magnetic field. We present details of the deduction of the 
time evolution of the quantum state of two nuclear spins 
manipulated by an external transverse magnetic pulse 
with some selective resonance frequencies, and show how 
to achieve an arbitrary non-adiabatic geometric quantum 
phase gate. Then, as a results, the controlled geometric 
phase gate as well as A-A phase are easily obtained. 

We first recall the definition of the quantum phase 
gate (QPG). Denote T as an unitary transformation 
evolving the input quantum state |^(0)) to output one 
|*(t)) = T|*(0)), where T is a 2" x2" matrix (n the num- 
ber of qubits). The T defines a quantum gate. The gate 



is called as a phase gate if Ty 
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where {fa : i — 1,2,....} are the phases that are de- 
pendent on the evolution path in parameter space as 
well as the dynamics. Thereby each phase fa includes 
both the geometric and dynamical part in general. But, 
the dynamical part of fa can be vanished in some condi- 
tion. In the case, the phase becomes pure geometric and 
the transformation T is called as a geometric phase gate 



2 



(GPG). If the evolution of the state is instantaneous, T 
the non-adiabatic GPG. In our present work, we study 
the implementation of non-adiabatic GPG. 

The Hamiltonian of two nuclear spins with a weak 
Heisenberg type interaction in a constant longitudinal 
magnetic field along z direction is 



Hi 2 ) 



xy 



jj(2) + #(2) 



(wicrf + UJ2CJ2 + Jvl^) 
{Jafaf + Jal4), 



(1) 
(2) 



where isotropic coupling is assumed, lo\ and u> 2 are the 
Larmor frequencies of two nuclear spins, J the coupling 
constant, {erf, of, erf : i = 1, 2} the Pauli matrices, and 
h — 1. In the experiments, two different nuclear spins 
are selected, cl>i L02, we assume cji > o->2 and the longi- 
tudinal constant magnetic field is in the order of ITHz, 
so u>i,u>2 are much large than J and r\ — J /{lo\ — L02) << 

1. The experiment numbers of these parameters will 

(2) 

be given in late. Hx y is non-diagonal in a z represen- 
tation and gives the quantum fluctuation which yields a 

correction of order r] 2 , hence it can be ignored. Thus, 

(2) 

the Ising part of the Hamiltonian is a well pre- 

cise approximation. ij( 2 )(~ H^) has four eigenstates: 
{|00), |01>, 1 10) , |11)}, where denotes the spin up and 1 
the spin down. These states correspond to following 4 
eigenvalues respectively: 

ex = — (wi + u>2 + J)/2, £2 = — (tjJi — — J)/2, 
£3 = (^1 - ^2 + J)/2,e 4 = {ll>i + uj 2 - J)/2 

In this paper, we want to find out an appropriate mag- 
netic pulse to achieve arbitrary non-adiabatic geometric 
QPG. Therefor a rectangular transversal magnetic pulse 
with 4- frequency (fii, f^, ^3, ^4) is applied to the sys- 
tem. The Hamiltonian becomes 



where 



Htotal = + Hpulse(t), (3) 



H pulse (t) = -lj2^h ik (t), (4) 
i=l fe=l 

h ik (t) = h ik (f k (t)at + f k {t)cjr). (5) 



H pu i se (t) is the external magnetic pulsed where f k (t) = 
e i(n fc t+*fc) ) a ± — a x ±i a y ^ (jf ailc l a V are Pauli operators. 

h\ k — J\h k and fi2 k — 72/i* and /i£ the amplitudes of 
magnetic pulse. 71 and 72 are the gyromagnetic ratio 
of two different nuclear spins. It is assumed that the 
magnetic pulse is applied in the duration [0, r] : 

h k ^ 0,<e [0,r] 

h k = 0,t$ [0,r],k = 1,2,3,4. 



a 



n> 



10: 



01: 



00> 



FIG. 1: The sketch of the energy levels of the Hamiltonian 
Htotal and two multi-photon transition schemes 



The phases {$&} and duration r of the magnetic pulse 
will be designed below in detail. Set 



n 3 



£3 - ei 
£4 - £3 
e 2 - ei 



: LUX 

: UJ 2 



J, 
J, 
J, 



SI4 = £4 — e 2 = uii — J, 



|00) 
|10) 
|00) 
|01) 



Ho) 

111), 

101), 

111). 



(6) 
(7) 
(8) 
(9) 
The 



Corresponding transitions are shown in figure. 
Schrodinger equation is 

ihd\v{t)))/dt = H tota i(t)\v(t)), 

The solution can be expended as follows 

!*(*)> = I>(t)K) (J = 1,2,3,4) (10) 

3 

where |mi), \rri2), ITO3), (^4) respectively corresponds to 
the state |00), |01), 1 10) , |11). When the system is initially 
at the state |^(0)) = 1 00) , the pulse with 4 frequencies 
{tti : i — 1,2,3,4} leads a coherent evolution of four 
states |00),|01),|10),|11). Since those transitions satisfy- 
ing resonant condition are dominating, we can ignore all 
other minor detuned transitions. This is so-called reso- 
nant transition approximation (RTA). In RTA, the states 
transferred each other by the magnetic pulse are called 
as the active states, the rest the inactive. When we apply 
an magnetic pulse with 4 frequencies (see Fig.l), in RTA 
the Hamiltonian H to tai(t) can be well approximated by 
following H(t), 



H(t) = 



V 



■^/l*(t) 





-¥/ 3 (t) -^h(t) 



£2 




£3 



4/a(*) 



ffm -fm) £ 4 j 



(ii) 

where h\ = hu, h 2 = h 2 2, h 3 = hi 3 , h± = /i 24 and 
is the complex conjugate of f k (t). 

Now we want to find out the exact solution of Hamilto- 
nian without adiabatic approximation. We do an unitary 
transformation U(t) as follows 



(12) 
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From Schrodinger equation ihd(\^(t)))/dt = H{t)\V{t)}, 
we obtain an equation of \ip{t)) : 

ihdmt)))/dt=H rot (t)m)), (13) 

where 

H ro t(t) = U{t)H{t)U + {t) - iU{t)dU + (t)/dt. (14) 

It is assumed that the matrix U{t) is diagonal: Uij(t) = 
5ije % (' Pit+0i ') where {<Pi,6i} are some parameters to be 
determined. When the parameters {<fii,6i} satisfy the fol- 
lowing relations: 



Ql 


= 


- <Pl,®l 


= e 3 -0i 


(15) 




= ¥4 


- f3, $2 


= 64 — 3 , 


(16) 


n 3 


= ¥>2 


- ip u $ 3 


= 02—01 


(17) 


n 4 


= <P4 


- lf 2 , $4 


= 04 — 02 


(18) 



It can be easily found that 
H rot (t) = 



2 





V 



-<pi - }J f 

\ e 2 - V2 

\ e 3 -cp 3 - 

£4 



\ 



2 

h2 

2 



(19) 



which is time independent. We simply denote H rot (t) as 
H rot . Due to Oi = 63 - ei, O2 = £4 - £3, ^3 = £2 - ei, 
CI4 = £4 — €2, then considering the relations of (15)-(18), 
we get 



£2-^2 = ^3 + £1 - <^2 
£3-^3 = fii + ei - ip 3 

£4-^4 = ^3 + ^4 + £l 



£l - <Pl, 
£l - <Pl, 

if4 = ei - pi- 



Setting tpi = £1, equation (19) becomes 

H ro t(t) = H rot 
From equations (15)-(18), we have 








h 3 
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hi 
2 




h 3 
2 










hi 
2 








I 





hi 
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h 2 
2 



hi 

2 

_h2 

2 

J 



(20) 



lfi2 = Q 3 + (fii,lfi 3 = (fii, 

lf4 = Ol + O2 — <Pl, 

02 = $3 + $1,03 = $1+01, 

4 = $l+$2-01- 

Four eigenvalues of H rot can be found as following 

e 1 = ^^Ja + VWT^, 



E2 
E 3 
E 4 



4 
V2 
4 

-V2 



A + VB + C, 

= -I^a-VbTc, 



where 



A = ^l B = J2ht, 

»=i i=i 

C = 8hih 2 h 3 h 4 + 2/i 2 /j| + 2/j 2 /iJ 

we further set hi = hi and /i 2 = h 3 . It results 

1 1 

E x = --{hi +h 2 ), E 2 = -{hi +/12), 

S3 = -^(/ll-te), ^4 = ^(/il - /l 2 ). 



(21) 
(22) 



Four corresponding cigcnfunctions are = 
^2j = iCij\m,j),i = 1,2,3,4. The coefficients matrix 



C and its inverse C 1 are 

C = C~ l = 



( l 1 1 1 \ 

1-1-1 1 

1-11-1 

Vi 1 -1-1/ 



(23) 



\ m i) = Ej=iCy = 1,2,3,4.. If the system is 

initially at the one of the eigenstates: |W^)(0)) = \rrii), 
{i = 1,2,3,4), the state of system at time t will be 

|*«(t)) = U+{t)exp[-iH rot t}\ mi ) 
4 

4 

= S C^Vfc)WI^ + (*)l^)exp(-i^i) 
4 

= E ^V^WexpH^C^K) 
= 24 i) (t)|m fc ),t = 1,2,3,4, (24) 



where 



4°w 



2 exp(-i£7 J -t)C J - fc ) e -*(^'+«») 
4 <) (t)e- i(vfe * +0fe) (25) 
£C^< "•'''>• (26) 



The initial condition is |cf } (0)| 2 = 1 and |c£° (0)| 2 = 
{k + *)■ 

For a cyclic evolution after time r, 



|* W (r)) =e tft M|*(0)), 



(27) 



where the total phase is /3,(t) = arg(*W(0)|*W(r)) 



S${t) + Sq'{t). Sq'{t) is the geometric phase. The 
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dynamical phase 6^ (r) can be calculated by formula: 
5®(t) = I dt{^(t)\H(t)\^(t)) 



dt^ C W\t)c):>(t)(m k \H(t)\rn k ,l28) 

k,k' 



In appendix A, we obtain 



= ^(t)=4 8) (*)=^ 4) (*) 



'h \ (hi 

cos [ 1 COS 1 

2 I 2 



= 2f(t)=ci 8) (*)=#>(*) 
= iC os[yt]sm[yi] 

= 4 2) w = ^i 3) w = 4 4) w 

= *sin[yt]cos[yt] 

. r^l.i . r ^2,i 

= sm [y*] sm[y t] 



(29) 
(30) 

(31) 
(32) 



When t = t and satisfies h±T = 2to7t, h^r = 2nn, namely 

n 



t = 2mir/hi, hi 
we have 



in 



-h\., m,n — 1, 2, 3, . 



(33) 



1 :i = 1,2,3,4, 
: i± j. 



It directly yields 

= Y. e ' i{Vht+9h) '^\ t )\ m k) 



= e- i ^* T+e ^\m i ),i = 1,2,3,4. 



Equation (33) is the condition to achieve a phase gate. 
In appendix B, we give the proof of that the dynamical 
phases also vanish at same condition (33): 



5«(r) = 0,i=l,2,3,4. 



(34) 



Hence, for any initial state |^(0)) = X^ c *(0)l m i) l ne 
wave function |*(r)) is 



-ie. 



|*(t) = A 



e 

\ 





ie 2 



\ 






l*(0)>- (35) 



Where A(= cos(m7r/2) cos(n7r/2), A — ±1) only con- 
tributes a global phase. The phases {Qi = ipir + : 
i = 1,2,3,4} are all geometric. Hence a non-adiabatic 
geometric phase gate is achieved. 



From equations (15) to (18), there are relations ip 2 — 

^3 + <Pl, ¥>3 = Ox + ^X, </?4 = fix + ^2 - </>X,</>X = e X- 

It is similar to have # 4 = $i + $ 2 - #Xj #2 = $3 + $x, 
#3 = $1 + #1. Since all {fii} and e\ are fixed, all {<^i} 
are changeless. are the phases of four harmonic 

waves of the external magnetic pulse and can be freely 
adjusted. Hence we can use four free parameters to 
change the values of {0i : i — 1, 2, 3, 4}, then to alter four 
phases Oj(= </?i r + 0») m gate individually. As a result 
the arbitrary non-adiabatic geometric phase gate of two 
qubits is implemented by means of NMR. The duration 
time r of magnetic pulse is determined by equations (33): 
t = 2rmr/hi,h2 — —hi where hi and h 2 are related 
to the amplitude of magnetic pulse and are adjustable. 
Taking m = n = l,it gives hi = h 2 and r = 2n/hi. 

Set 61 — —(fir and 8 2 = —^2^ in equation (35), we ob- 
tain an arbitrary geometric controlled phase gate (CPG) 
by changing the parameters {#3, #4}: 



T, 



CPG 



A 



/ 1 
1 

e~ i( - V3T+03 ^ 



\Q 











-i(ip4T + 8 4 ) 



\ 



, A = ±1. 



Further, if we adjust the parameters {6i} to satisfy fol- 
lowing conditions 



IflT ■ 



2,3,4. 



The wave function at time r is recovered to initial state 
with a phase difference only. 



|*(r))=±e" 



l*(o)>. 



Where ipiT + 81 is just the A- A phased. 

Finally, we give brief discussion on the validity of Ising 
model and our RTA and some estimation of the time r to 
achieve the geometric arbitrary phase gate. It is known 
that the nuclei 1 i7, 13 C, 15 A, 19 F, 31 P,etc. have spin S = 
1/2. Taking the 1 H and 13 C as an exampl e - 27 ' 29 , the Lar- 
mor procession frequency of nucleus 1 H : oji ~ 500 MHz 
at field Bq = 1.87T. The gyromagnetic ratio of 13 C has 
only 25% of 1 H U , so u 2 ~ 125MHz. uj y -^ 2 ~ 375MHz. 
Assume the coupling J ~ 200iJ g 27 i 29 (in the references 
J - 200Hz), we have J/{u x - uj 2 ) ~ 0.54 x 10~ 6 (<< 1). 

(2) 

Hence the quantum fluctuation from H X y can be ignored. 
It shows that our Ising model Hamiltonian can well de- 
scribe the physics. 

We set the intensity of the external transverse mag- 
netic pulse hi ~ 10~ 2 T, it yields h x ~ 2.8MHz. From 
the equation (33), we have 

r = 2rnn/hi ~2.2xrax 10 -6 sec . 

Meanwhile, h 2 should be adjusted to meet the condition 
h 2 = —hi. Since gyromagnetic ratio of 13 C is 0.25 of 

the one of 1 H, the amplitude h 2 must be four times of hi 
(h 2 ~ 4hi), then the condition h 2 = zjrhi can be satisfied. 
We choose the one of solutions, n = 1, m = 1, that yields 
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h 2 {= hi) ~ 2.8MHz and r = 2tt//ii - 2.2 x 10~ 6 sec. 
Due to f2i : 4 = wi ± J, f22,3 = 0J2 T >A we have fijT ~ 
2.5 x 10 2 — 1.0 x 10 3 >> 1 for all i. Larger m or smaller 
h\ yields more longer r. Thus, the precision of our RTA 
can be well ensured. 

The phase decoherence time of nuclear spin is the 
order of 10 4 seo 3 -, but the operating time r of such 
non-adiabatic geometric phase gate is only the order of 
10~ 6 sec . The ratio r^/r(~ 10 10 - 10 9 ) in NMR is high. 

To summarize, We have studied Ising model of two 
nuclear spins. A single tranverse rf magnetic pulse with 
multi-frequency is applied to manipulate quantum state 
simultaneously. By means of RTA, an arbitrary non- 
adiabatic geometric phase gate. As a result, the con- 
trolled phase gate as well as A-A phase are also ad- 
dressed. From estimation of parameters, we believe that 
our RTA is a good approximation and the operation time 
of GPG t is order of 10~ 6 — 10~ 5 sec. The absence of 
dynamical contributions to the phase gate can decrease 
the error of dynamical fluctuation in quantum informa- 
tion process. Meanwhile, the non-adiabatic manipula- 
tion of quantum state in practice quantum computing is 
required that leads the study of non-adiabatic manipula- 
tion to implement the geometric phase gate. 
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helpful discussion. This work is supported by the China 
National Natural Science Foundation (Nos. 10234010 and 
10674027) and by 973 project (No.2002CB613504) of 
China Ministry of Science and Technology. 



Appendix A 



In this appendix, we would prove |c^(r)| = 8\ 
when h\T — 2rmr and h 2 r = 2nn where m and n 
are integer numbers. According to eq.(26), cj|. (i) = 

4 

J2 <",,'< •'<"„,. and eq.(23) of matrix C, we get 

3=1 



3=1 

— _|g— i-Eit , e ~iE 2 t e —iE 3 t , e ~iE 4 t 

4 

If .hi + h 2 . M-hi 

= 2{ COS{ ^2— )+COS{ ^2— ) 

'hi \ (h 2 

cos 1 COS 1 

2 V 2 



} 



e — iEit _ e -iE 2 t _ e ~iE 3 t _j_ e ~iE 4 t-i 



3=1 
1 

4 

i (hi + h 2 )t . (hi - h 2 )t 

-sin sin 

2 L 2 2 J 

hi h 2 
% cos[— 1\ sin[— 1| 



It is similar to obtain following results: 

= c?> (*) = <£»> (*)=#>(*) 



l COS 1\ Sill 1\ 

1 2 1 1 2 1 



= ism[yt]cos[yi] 

#>(t) = if (t)=4»> (*)=#>(() 

. r hi , . r h 2 
= sm[y t] sm[y t\ 

When t = t, hiT = 2miT and h 2 T = 2mr (m, n : integer 
numbers), following results are obtained: 

\^\r)\=l;i j \r) = 0,i^j. 
Appendix B 
In this appendix, we want to prove the equation (28): 



*(<) 

J D 



0:i= 1,2,3,4. 



When hi = h±, h 2 = /13, the definition of H(t) becomes 
(see eq.(ll)) 



H(t) = 



( ex -^/ 3 (*) -^h{t) \ 

-¥/s*(*) e 2 -^/ 4 (t) 

-^m) e 3 -^/a(t) 

v -*fft(t) -ifm) e 4 , / 



where fk(t) = e l ( f2fc * + *'" . According to the definition of 
dynamic phase, eq.(28), we have 



c(0 



dt^cW*(t)c%(t){m k \H(t)\m k , 



k,k' 



J ° k=l 



- |zEg< i )*(^2( i )^) e -»(^3t+*3)-«(vi*+ei)+i(v2t+6'2) 



+ (t)cf' (£} e -i{Clit+'bi)-i < ,V2t+e 2 )+i( Vi t+8 i ) 

+ _c|W* (t)($ (t) e - i ( n2t +* 2 )~ i (^ 3t+03 )+ i (^ 4t + 04 ) 
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Using the relations of eq.(15) — (18), above equation be- 
comes 

= [ dtjy?\t)^{t)e k 
Ja fc=i 



+ y pf*^ (*) + Q?>* (t)^° (*)] + fc-c.}. (Bl) 



We do following calculations: 
Jo 



^- / + cos(/i 3 t)][l + cos(M)] 
r + y (it cos(h 3 t) cos(/i 4 i)^ 

y dt(cos(ti3 + hi)t + cos(/i 3 — /i4)i 



4 



Input ei = — (wi + u)2 + J)/2, £2 = — (wi - W2 — J)/2, 
£3 = (a>i — CJ2 + </)/2, and e 4 = (wi + o>2 — *0/2 to above 
equations, we get 



f dtjy^\t)^\t)e k = T - = 0. (B2) 

Ja fe=i \fe=i / 



Moreover 



2 



Jo 

— / dt{—i sinf— t] cos[— t] cos[— t] cos[— t] 
2 Jo 1 L 2 J L 2 J L 2 J L 2 J 



+ i sin[y i] sin[y i] sin[y i] cos[y i]} 



/11 r 



dt{—ism[hit] cos[/i2*]} = 



(B3) 



dt^ 1) *(t)^ 1) (t)e 2 



= ^ / dt[l + cos(/iit)][l-cos(ft 2 t)] 
4 Jo 

y dt cos(hit) cos(h2t)^j — ~j~ T - 



£2 
4 ' 



It is similar to have 



/' 

= e 3 y disin 2 



2 / "2 



COS I —t \ = -T, 



£3 



/ 

= £4 / 

Jo 



(// <in I ytj sin fytj = jT. 



2 



Jo 

= y y s in[y t] sin[y t] sin[y t] cos[y t] 
- * cos [y *] sin[y t] cos[y i] cos[y t]} 

1, rr 

= -4 / dt{isin(/i 2 t)cos(/nt)} = (B4) 
2 Jo 



Inserting results of (B2-B4) to (Bl), we obtain S^\t) = 

(2, 3 4) 

0. Similar calculations also can 

show<5)y J ^(T) =0. The 
conclusion of S^\t) = is proved for all i. 



1 P. W. Shor, In proceeding of the 35th Annual Symposium 3 David P. DiVincenzo, Phys. Rev. A 51, 1015 (1995). 

on the Foundations of Computer Science (IEEE Computer 4 M. A. Nielson and I. L. Chuang, Quantum Computation 

Society, Los Alamitos, CA 1994) P124. and Quantum Information Cambridge University Press, 

2 Lov K. Grover, Phys. Rev. Lett. 79, 4709 (1997). Cambridge, UK, 2001. 



7 



5 Seth Lloyd, Phys. Rev. Lett. 75, 346 (1995). 

6 Adriano Barenco, Charles H. Bennett, Richard Cleve, 
David P. DiVincenzo, Norman Margolus, Peter Shor, Ty- 
cho Sleator, John A. Smolin, and Harald Weinfurter, Phys. 
Rev. A 52, 3457 (1995). 

7 Seth Lloyd, Science 261, 1569 (1993). 

8 Q. A. Turchette, C. J. Hood, W. Lange, H. Mabuchi, and 
H. J. Kimble, Phys. Rev. Lett. 75, 4710 (1995). 

9 D. Deutsch, A. Barenco, A. Ekert. Proceedings of the 
Royal Society of London Series A-Mathematical and Phys- 
ical Sciences 449 (1937). 

10 M. V. Berry, Proc. R. Soc. London, Ser. A 45,392 (1984). 

11 A. Shapere and F. Wilczek, Geometric Phases in Physics, 
(World Scientific, Singapore, 1989). 

12 Frank Wilczek and A. Zee, Phys. Rev. Lett. 52, 
2111(1999). 

13 Ferdinand Schmidt-Kaler, Hartmut Haifner, Mark Riebe, 
Stephan Guide, Gavin P. T. Lancaster, Thomas Deuschle, 
Christoph Becher, Christian F. Roos, Jiirgen Eschner & 
Rainer Blatt, Nature (London) 422, 408 (2003). 

14 D. Leibfried, B. DeMarco, V. Meyer, D. Lucas, M. Bar- 
rett, J. Britton, W. M. Rano, B. Jelenkovic, C. Langer, T. 
Rosenband & D. J. Wineland, Nature (London) 422, 412 
(2003). 

15 A. Rauschenbeutel, G. Nogues, S. Osnaghi, P. Bertet, M. 
Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. 
83, 5166 (1999). 

16 Shi-Biao Zheng, Phys. Rev. A 70, 052320 (2004). 

17 R. Garcfa-Maraver, R. Corbalan, K. Eckert, S. Rebic, M. 
Artoni, and J. Mompart, Phys. Rev. A 70, 062324 (2004). 

18 Jonathan A. Jones, Vlatko Vedral, Artur Ekert and 
Giuseppe Castagnoli, Nature (London) 403, 869 (2000) 

19 Ralf Stadelhofer, Dieter Suter and Wolfgang Banzhaf, 
Phys. Rev. A 71, 032345 (2005) 



20 Bin Zhou, Ruibao Tao and Shun-Qing Shen, Phys. Rev. A 
70, 022311 (2004) 

21 Yu Shi and Yong-Shi Wu, Phys. Rev. A 69, 024301 (2004). 

22 Y. Aharonov and J. Anandan, Phys. Rev. Lett. 58, 1593 
(1987). 

23 Go Yusa, Koji Muraki, Kei Takashina, Katsushi 
Hashimoto and Yoshiro Hirayama, Nature (London) 434, 
1001 (2005) 

24 Isaac L. Chuang, Lieven M. K. Vandersypen, Xinlan Zhou, 
Debbie W. Leung and Seth Lloyd, Nature (London) 393, 

^ 143 (1998). 

25 Lieven M. K. Vandersypen, Matthias Steffen, Gregory 
Breyta, Costantino S. Yannoni, Mark H. Sherwood and 
Isaac L. Chuang, Nature (London) 414, 883-887 (2001). 

26 Jonathan A. Jones, Michele Mosca and Rasmus H. Hansen, 
Nature (London) 393, 344 (1998). 

27 M. A. Nielsen, E. Knill, and R. Laflamme. Nature (Lon- 
dond) 396, 52 (1998). 

28 S. Somaroo, C. H. Tseng, T. F. Havel, R. Laflamme, and 
D. G. Cory, Phys. Rev. Lett. 82, 5381 (1999). 

29 E. Knill, R. Laflamme, R. Martinez, and C.-H. Tseng. Na- 
ture (London) 404, 368 (2000) . 

30 L. M. K. Vandersypen, Experimental Quantum Computa- 
tion with Nuclear Spins in Liquid Solution, Ph.D. thesis 
Stanford Universityd (2001). 

31 Wurong Zhang and D. G. Cory, Phys. Rev. Lett. 80, 1324 
(1998). 

32 Garett M. Leskowitz, Nima Ghaderi, Ryan A. Olsen, and 
Leonard J. Mueller, J. Chem. Phys. 119, 1643 (2003). 

33 Michael N. Leuenberger, Daniel Loss, M. Poggio and D. 
D. Awschalom, Phys. Rev. Lett. 89, 207601 (2002) 

34 R. D. Inglis, Phys. Rev. 58, 577 (1940) 



